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Abstract: [ Aim] Anopheles stephensi is a major urban human malaria vector in South-East Asia 
contributing about twelve percent of malaria cases in India. The gene flow in An. stephensi belonging to 
the areas under Enhanced Malaria Control Project (EMCP) and non-EMCP of north-eastern regions of 
Madhya Pradesh (central India) was studied as initially there was a decline in malaria cases in EMCP 
areas due to various malaria control strategies but malaria appeared to be back-track with malaria cases 
increasing readily, indicating that the overall malaria risk remained stable. [Methods] The An. stephensi 
mosquitoes collected from four different areas under EMCP and non-EMCP from north eastern Madhya 
Pradesh ( central India) were genotyped using 7 microsatellite loci for analysis of various population 
parameters. [Results] The markers were found to be highly polymorphic in all the populations. Not 
much genetic diversity was found between the two regions. Higher genetic differentiation was observed in 
eastern populations under EMCP ( Fẹ = 0.0485, Rẹ =0. 1112) than northern populations under non- 
EMCP (Fo, =0.020, Rey =0.0145) and high gene flow (12.90, 6.16, 5.06 and 2.38) was observed 
between EMCP and non-EMCP areas. The higher sensitivity of Ry, than Fr indicated that the 
differentiation was probably caused by the mutation but not due to genetic drift. [Conclusion] The study 
revealed substantially high gene flow within and between EMCP and non-EMCP areas. The high gene flow 
combined with development of insecticide resistance seems to be a potent reason for increase in malaria 
cases in EMCP as well as non-EMCP areas. 
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1 INTRODUCTION 


Malaria continues to be a major public health 
India 
substantially to global malaria incidence ( Hay et al., 


problem worldwide and contributes 


2010). Despite more than five decades of intensive 
control efforts, malaria is still a challenge to the 
Indian health system. Anopheles stephensi is an 
important urban malaria vector in Indian 
subcontinent contributing about twelve percent of 
malaria cases in India ( Adak et al., 2005) and its 
has 


resistance. 


control been hampered by insecticidal 
There is no information available about 
the population structure and gene flow of An. 


stephensi in India, except in Rajasthan ( Vipin et al., 





2010a) and in Haryana ( Vipin et al., 2010b). The 
population genetics studies on An. stephensi using 
different molecular markers in different parts of the 
world have been reviewed by Gakhar et al. (2013). 

Microsatellites are highly polymorphic molecular 
markers most widely used in population genetic 
studies because of their high degree of variability due 
to high mutation rate. Moreover, these are 
selectively neutral and variations are found even in 
populations showing low level of allozyme and 
mitochondrial variation. 

Madhya Pradesh ( MP ) is the largest state 
situated in the central India (308 252 km’) having 
50 districts with forest cover of about 31% with 
intense perennial malaria transmission in these areas 


(Mishra et al., 2012). 


Malaria is complex in 
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Madhya Pradesh because of vast tracts of forest with 
tribal settlement (Singh et al., 2004). Enhanced 
Malaria Control Project (EMCP) was introduced in 
1997 by National Anti-Malaria Programme (NAMP) 
with World Bank assistance to control the malaria in 
tribal/endemic areas. Initially, there was decline in 
malaria cases but it could not be maintained and it 
appeared to be 
increasing readily , indicating that the overall malaria 
risk remains stable (Singh et al., 2009). In the 
present study, we have therefore studied the genetic 
structure of An. stephensi in the areas under 
Enhanced Malaria Control Project (EMCP) and non- 
EMCP of central India using microsatellite DNA 
markers developed by Verardi et al. (2002). The 
information on An. stephensi between these regions 
should help study epidemiology. The knowledge of 
population genetic 
implementation of 


back-track with malaria cases 


structure is prerequisite for 


adequate malaria control 
strategies, particularly the transgenic technology. 
This study should also help to understand the reasons 
for increase in malaria cases, particularly in highly 


endemic areas, where EMCP was implemented. 


2 MATERIALS AND METHODS 


2.1 Sample collection and identification 

The adult female An. stephensi mosquitoes were 
collected from non-Endemic/non-EMCP, i. e., 
Gwalior (26°12'46. 54”N, 78°11'0. 66"E), Sagar 
(23°49'51.66"N, 78°44’41.16”E) , and Endemic/ 
EMCP areas, i. e., Umaria (23°31'31. 13”N, 80° 
50’28.18”E) and Anuppur (23° 6'5.95”"N, 81°41’ 
52.49” E ), in north-eastern region of Madhya 
Pradesh (Fig. 1). 

Mosquitoes were 


hand 


indoor 


collected using a 
held suction tube and torchlight from 
resting sites. The samples were collected from two or 
three different areas around 10 km apart. Distance 
between different collection sites is given in Table 1. 


Intact mosquitoes were morphologically identified 





Fig. 1 Map showing collection sites of 
Anopheles stephensi mosquitoes 
Star indicates non-EMCP areas and dotted circle indicates EMCP areas. 


Green circle represents Vindhya mountain range. 


(Nagpal et al., 2003) and were preserved in 95% 
ethanol and stored at —20°. 
Table 1 Geographic distances between population 
pairs of Anopheles stephensi 


Population pairs Distance (km) 





Gwalior-Sagar 270 
Gwalior-Umaria 400 
Gwalior-Anuppur 495 
Sagar-Umaria 216 
Sagar-Anuppur 310 
Umaria-Anuppur 95 


2.2 Genomic DNA extraction and microsatellite 
genotyping of samples 

Genomic DNA of 30 mosquitoes was extracted 
from the legs and wings of individual mosquitoes and 
stored at — 20° (Vipin et al., 2010a, 2012b). 
Seven microsatellite loci developed by Verardi et al. 
(2002 ) were used for genotyping 30 individuals 
from each population of An. stephensi. Each pair 
of microsatellite markers were labeled by 6-FAM. 
Multiplex PCR was by grouping 
the microsatellite according to their 


carried out 
primers 
allele size range and annealing temperature ( Table 
2). The reaction conditions were initial heat 
activation at 95°C for 15 min, and denaturation 


at 94°C. 


Table 2 List of microsatellite loci and their combinations used for multiplex PCR 





Locus Forward sequence (5' —3’) Reverse sequence (5’—3") Combinations used T, 
D8T TCCTTCACAGCTCGGAACTT AGAGTATCGGCTTATGGGATG A7 & A10 60C 
A7 CACAAATCAACCCCCTTCAT CTTCAATCACCGAAACAACAA H2II & E12 59% 
A10 GCTGGACACGCAAGTAGG ATAAAAGGGCACGAAACTCC F10 & B2N 59% 
H211 ACAATCGGTCGGAGTTCATC AAAGCAGAGCCAAATCAAAAG D8T 58C 
E12 ACATGCAACACCGAAACACC CCACAGGAGGACGACTTGG 
F10 GTGGGTGCAGTAGTGGGATT GATCCAAACCGTGGTACAGC 
B2N GCTGGGAAGTGCTACAGAGG AAGGTTCGAGTACGGACGTG 


T, = Annealing temperature for 30 s, annealing at 57 -60°C depending on combination used and extension at 72°C for 90 s. The numbers of cycles used 


were 30 —45. Genotyping was done by ABI 3730xl Genetic Analyzer and the data were analyzed by the Genemapper 4. 0 software. 
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2.3 Data analysis 
Genetic populations was 


measured by number of alleles, 


Diversity within 
observed and 
expected heterozygosity and effective number of 
alleles at each locus using POPGENE version 1. 32 
(Yeh et al., 1999) and GeneAIEx version 6. 41. 
Fisher exact test available on Genepop version 3. 1 
(Raymond and Rousset, 1995) was used to test the 
Hardy-Weinberg Equilibrium in each population. 
Exact test of genotypic disequilibrium between all the 
pairs of loci in all the populations was performed to 
know the independence of loci using Genepop 
version 3. 1 (Raymond and Rousset, 1995). Micro- 
checker computer programme ( Oosterhout et al., 
2004) was used to estimate the frequency of null 
alleles in all the populations at each locus. The 
population bottleneck was studied by Bottleneck 
version 1. 2. 02. The 
populations was measured by F-statistics and R- 
statistics by FSTAT version 2. 9. 3. 2 ( Goudet, 
1995 ). 


genetic differentiation and geographic distance, the 


differentiation between 


To investigate the relationship between 


regression of F (1 - Fr) on the natural log of 
geographical distance was studied (Fo, 1997). The 
dendrogram was constructed using the Nei’ s genetic 
distances ( Nei, 1978) by UPGMA method in 
POPGENE computer programme. 


3 RESULTS 


3.1 Genetic of An. 
populations 


Genotypes of four populations ( Gwalior, Sagar, 


variability stephensi 


Umaria and Anuppur) scored at 7 microsatellite loci 
indicated that these were highly polymorphic in all 
the populations with number of alleles ranged from 5 
to 19 (Table 3). Effective number of alleles ranged 
from 2. 00 to 12. 24. A total of 87 alleles were 
observed in all the populations. The E12 locus was 
found to be the most polymorphic with 22 alleles. 
This locus had a maximum of 19 and 16 alleles in 
northern populations (non-EMCP) , while 13 and 11 
The observed 
heterozygosity varied from 0. 333 to 1. 000 among 


in eastern populations ( EMCP ). 


seven loci and the expected heterozygosity varied 
from 0.56 to 0.93 (Table 3). The number of alleles 
(N, ) 
appeared to be slightly higher in the northern 
populations under non-EMCP (10.21, 0.77) than 
eastern population (8.4, 0. 75) under EMCP of 
Madhya Pradesh. The four populations ( Gwalior, 


and the expected heterozygosity ( H, ) 


Sagar, Umaria and Anuppur) showed more or less 
similar allelic distribution pattern (Fig. 2). 


Table 3 Genetic diversity at 7 microsatellite loci in 
Anopheles stephensi populations 



































Populations 
Locus 
Northern (non-EMCP ) Eastern ( EMCP) 
Gwalior Sagar Umaria Anuppur 
(N =30) (N=14) (N =30) (N =24) 
N, 9 6 5 7 
N, 5.21 3.18 2.75 3.21 
D8T H, 0. 733 0.566 0.566 0.625 
H, 0.82 0.69 0. 64 0.70 
Fis 0. 109 * 0. 189 * 0. 126 * 0.114* 
N, 11 10 9 9 
N, 4.42 4.82 5.50 5.16 
A7 H, 0.766 0.666 1. 000 0.625 
H, 0.78 0.80 0. 83 0. 82 
Fig 0. 026 0.177* -0.205 0.245 * 
N, 5 5 7 5 
N. 2.83 3.04 2.76 3.51 
A10 H, 0. 600 0.615 0. 800 0.583 
H, 0.65 0.68 0.65 0.73 
Fis 0. 090 0. 103 -0.236 0. 204 * 
N, 11 8 10 9 
N, 4.13 2.22 5.37 3.72 
H211 H, 0. 633 0. 433 0. 333 0.583 
H, 0.77 0.56 0. 82 0.74 
Fis 0.181 0.229 0. 601 * 0.223 * 
N, 16 19 13 11 
N. 12.24 12.24 7.14 5.78 
E12 H, 0.700 0. 600 0. 500 0.708 
H, 0.93 0.93 0. 87 0. 84 
Fis 0.253 * 0. 361 * 0. 432 * 0. 164 * 
N, 14 16 8 14 
N. 9.32 11.76 3.05 8.53 
F10 H, 0. 866 0.966 0. 466 0.875 
H, 0.90 0.93 0.68 0.90 
Fis 0.046 * -0.039 0.321 * 0.030 
N, 6 7 5 6 
N, 3.18 3.21 3.62 3.27 
B2N H, 0. 400 0. 433 0. 800 0. 666 
H, 0.69 0.70 0.73 0.70 
Fis 0.430* 0.385 * - 0. 087 0.061 
N, 10. 28 10. 14 8.14 8.71 
N, 5.91 0. 578 4.31 4.74 
Mean 
H, 0.671 0.611 0. 638 0. 666 
H, 0.79 0.75 0.74 0.77 
N, 1.5 2.017 1.079 1.169 
D N, 1.322 1.633 0. 644 0.735 
H, 0.056 0. 068 0.089 0.039 
H, 0.038 0.051 0.035 0.028 





Average N, =10.21, 
Average H, =0.77 


Average N, =8.4, 
Average H, =0.75 


* Significant at P <0.05; N = Number of samples from each population; 
N, = Number of different alleles; N, = Effective number of alleles [ Kimura 
and Crow, (1964) ]; H, = Observed heterozygosity; H 


e 


= Expected 
heterozygosity; Fıs = Inbreeding coefficient calculated according to Weir 


and Cockerham (1984). 
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Fig. 2 Allelic distribution pattern of 7 
microsatellites in Anopheles stephensi populations 
A = Number of different alleles; B = Number of different alleles with a 
frequency =5% ; C = Number of effective alleles; D = Number of unique 


alleles; E = Number of common alleles with a frequency <50% . 


Three alleles 88, 173 and 87 of D8T, A10 and 
B2N respectively, were found to be predominating 
populations, allele 87 
corresponding to B2N locus was found to be the most 
frequent and predominating ( Table 4). Some unique 
alleles corresponding to one or other populations 
were also observed but their frequency was very low 
(data not shown). More unique alleles were present 


in northern population ( non-EMCP) than the eastern 
population (EMCP). 


across all however , 


3. 2 Hardy-Weinberg equilibrium and linkage 
disequilibrium observed in An. stephensi populations 

The significant deviation from Hardy-Weinberg 
equilibrium (HWE) was found in 17 (P <0. 05) 
out of 28 tests within populations. Heterozygote 
deficiency was observed in all the cases of deviation 
from HWE. The loci D8T and E12 
significant heterozygote deficiency in all the four 
populations. The B2N, A7, H2II and F10 showed 


significant deficiency in two populations whereas 


showed 


A10 showed deficiency only in one population 
(Table 1). 

Significant linkage disequilibrium (LD) was 
observed only in 7 (8. 3%) loci pair out of 84 
pairs. These are D8T & A7, A7 &F10 and E12 & 
F10 of Sagar; D8T & A7, H2II & E12 of Umaria 
and A7 & H2II and D8T & B2N of Anuppur. When 
LD was studied the 
populations, only one locus pair ( D8T & A7) had 


locus wise-across all 
significant LD out of 21 pairs. 

3.3 Null alleles across microsatellite loci and 
population bottleneck in An. stephensi 


Null alleles the 
populations at 5 loci except at loci D8T and A10 


were present across all 
(Table 5). H „(expected heterozygosity under 
equilibrium ) based on two-phase mutation model 


(TPM) with fractions of multistep mutation set at 


Table 4 Alleles with maximum frequencies observed in Anopheles stephensi populations 


Alleles with maximum frequency (allelic frequency ) 

















i a D8T A7 A10 HII E12 F10 B2N 
Gwalior 98(0.283)  250(0.333) 173(0.517) 127(0.433) 188, 192, 194(0. 117) 178 (0. 200) 87 (0. 500) 
Sagar 88(0.433) 250(0.296) 179(0.481) 127(0.650) 194(0. 150) 171 (0. 133) 87 (0.450) 
Umaria 88(0.550) 252(0.300) 173(0.567) 141(0.267) 198 (0.267) 179 (0. 533) 85, 87 (0.350) 
Anuppur 88(0.479) 252(0.271) 173(0.438) 127(0.438) 194(0.271) 177, 179 (0. 188) 87 (0.435) 
Table 5 Null alleles across 7 microsatellite loci in Anopheles stephensi populations based on 
overall significant excess of homozygotes with their frequencies 
idi Populations 
Gwalior Sagar Umaria Anuppur 
D8T No (0.041/0. 048) No (0. 079/0. 094) No (0.052/0. 058) No (0. 038/0. 048) 
A7 No (0. 004/0. 004) No (0. 078/0. 086) No ( -0. 122/-0. 099) Yes (0. 112/0. 126) 
ALO No (0. 053/0. 038) No (0. 037/0. 043 ) No ( -0.197/-0. 112) No (0. 082/0. 101) 
H211 No (0. 099/0. 089) No (0. 119/0. 119) Yes (0.293/0. 418) No (0. 104/0. 112) 
E12 Yes (0. 115/0. 134) Yes (0. 171/0. 209) Yes (0. 207/0. 264 ) No (0.070/0. 077) 
F10 No (0.010/0. 014) No ( -0.030/-0. 027) Yes (0. 168/0. 180) No (0. 007/0. 004 ) 
B2N Yes (0. 192/0. 263 ) Yes (0. 191/0. 227) No ( -0.053/-0.049) No (0. 010/0. 020) 


The frequencies given were calculated according to the methods of Oosterhout et al. (2004) and Chakraborty et al. (1992). 


30% , 20% and 10% and also based on Stepwise 
Mutation Model (SMM) was calculated ( Comuet and 
Luikart, 1996). Significant values were found under 
SMM in population of Sagar, Umaria and Anuppur 


(Table 6). However, H, was less as compared to 


H 


eq’ 


distribution of alleles and a normal L-shaped curve 


There was also no mode shift in frequency 


was observed for all the four populations (Fig. 3). 
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Table 6 Two tailed probabilities for heterozygote excess 
and deficiency in Anopheles stephensi populations 
based on Wilcoxon signed rank test 


TPM with fractions of multistep 








Population mutations set to SMM 
30% 20% 10% 
Gwalior 0.8125 0.9375 0.9375 0. 10938 
Sagar 0.8125 0.9375 0.8125 0. 03906 * 
Jabalpur 0. 2968 0. 4687 0.4787 0. 2968 
Umaria 0.8125 1.000 0.9375 0. 03906 * 
Anuppur 0.8125 0.8125 0.8125 0.0156 * 


* Significant at P <0.05. TPM = Two phase mutation model; SMM = 


Stepwise mutation model. 


3.0 
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Allele frequency 
Fig. 3 Distribution of proportion of alleles in different allelic 
frequency classes in four populations of in 
Anopheles stephensi in Madhya Pradesh 


3. 4 Genetic differentiation and gene flow 
between populations of An. stephensi 

The level of genetic differentiation between 
pairs of populations was very high on eastern side, 
i. e., in the area under EMCP (F =0.0485). The 
average weighted Rover loci was 0. 111 (Table 7). 
The genetic differentiation between populations of 
northern region ( under non-EMCP ) was low as 
compared to eastern region. The average F',, value 
was 0. 020 and the weighted R,, was 0.0145. The 
gene flow was found to be higher in populations 
under non-EMCP (12. 25), i. e., Gwalior and 
Sagar, as compared to populations under EMCP 
(4.85), i.e., Umaria and Anuppur. 

The analysis of molecular variance (AMOVA ) 
of all the seven microsatellite loci confirmed 
the differentiation with the 
individuals, among individuals within populations, 


variation within 


among populations within groups and among groups 
being 78% , 17% , 3% and 2%, respectively ( Fig. 
4). AMOVA showed that the variation among 
populations within groups explained only 3% of the 
total variation while the variation from among 
individuals within populations and within individuals 
explained 17% and 78% of the total variation, 
respectively. 


Table 7 Genetic differentiation and gene flow ( N„) between populations of Anopheles stephensi 





Populations Gwalior Sagar Umaria Anuppur 
Gwalior 0. 000 0.0145 0.0879 0.0152 
Sagar 0.020 * (12.25) 0. 000 0.0777 0.0567 
Umaria 0.047 * (5.06) 0.095 * (2.38) 0.000 0.111 
Anuppur 0.019 * (12.90) 0.039 * (6. 16) 0.0485 * (4.85) 0.000 


* Significant at P <0.05 by bootstrapping over loci; F’gy values were calculated according to Weir and Cockerham, (1984) ; Rg; values were calculated 


according to Slatkin (1995) ; Nm =0.25(1 - Fs; )/F gy. Figures below diagonal are Fg, (N, 







EZA Among regions 

EZ Among populations 
ZZ] Among individuals 
S$ Within individuals 





Fig. 4 Analysis of molecular variance of microsatellite 
loci in Anopheles stephensi populations based on Fg, values 


3.5 Isolation by distance (IBD) in the An. 
stephensi populations 

The analysis of isolation by distance by 
performing Mantel test showed a significant and 


) values and above diagonal are Rg, values. 


m S 


negative (r= —-0.4618, P<0.05 ) correlation 
between F'g,/(1 — Fr) and the logarithm of distance 
(Fig. 5). 


4 DISCUSSIONS 


The genetic diversity between the populations of 
areas under EMCP and non-EMCP was slightly 
different as significant gene flow was observed 
between these areas. However, no previous study on 
population genetic structure of An. stephensi has 
been reported from Madhya Pradesh. The genetic 
diversity is measured with number of alleles and 
expected heterozygosity. The genetic diversity in 
northern region ( non-EMCP ) was slightly higher 
than the eastern population (EMCP). All the seven 
microsatellite markers used in the study were 
polymorphic with number of alleles at each locus 
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Geographic distance 
Fig. 5 Correlation between F,,/(1 — Fr) and logarithm 


of distance (in km) for pairwise comparison of populations of 
Anopheles stephensi at 7 microsatellite loci 


ranged from 5 to 19 (Table 1). The deviations from 
Hardy-Weinberg Equilibrium 
deficiency were observed in all the populations 


and heterozygote 


studied. The heterozygote deficiency results from the 
wahlund effect, inbreeding or presence of null alleles 
(Vipin et al., 2010a, 2010b). The wahlund effect 
or inbreeding might not be the cause of heterozygote 
deficiency as samples were collected from indoor 
resting sites only and mosquitoes belonged from a 
single region or two or three different sites within 5 to 
10 km range. Further, only 7 significant linkage 
equilibrium deviations were found ruling out the 
possibility of wahlund effect or inbreeding. The 
heterozygote deficiency seems to be caused by null 
alleles as confirmed by their presence at 5 loci across 
all populations. Null alleles, i. e., mutations may 
prevent annealing of primer resulting in either total 
failure of amplification or, most commonly 
amplification of single strand resting in scoring of 
genotype as homozygote (Lehmann et al., 1996). 
No evidence of mis-scoring of stutter was detected as 
there was no deficiency of heterozygote with alleles 
differing in size by a single repeat and an excess of 
large heterozygote. 

The significant deviation ( P<0.05 ) found 
under SMM and Anuppur 
populations (Table 2), appeared to be associated 
with H, less than H. indicating the absence of 
bottleneck ( Vipin et al., 2010a, 2010b). In a 


population that experienced bottleneck recently, the 


in Sagar, Umaria 


mutation-drift equilibrium transiently disrupts and 
the heterozygosity measured at majority of loci 
exceeds to HWE equilibrium heterozygosity (H.,, ) 
(Cornuet and Luikart, 1996). Normal L-shaped 


curve was observed as there was no mode shift in 
frequency of allele distribution, indicating that no 
recent event of population bottleneck has occurred in 
any population. 

High genetic differentiation was observed in 
Eastern populations (under EMCP ) than Northern 
populations (non-EMCP). High level of gene flow 
was also found between EMCP and non-EMCP areas 
except Sagar and Umaria, where stretches of 
Vindhya mountain range might be acting as a 
geographical barrier. The genetic differentiation was 
probably caused by mutation and not by genetic drift 
as indicated by the high sensitivity of Rg, than Fp. 
The isolation by distance was rejected by the 
negative correlation between genetic differentiation 
and geographic distance between sampling sites as 
evident by mantel’ s test. This indicated that 
neighboring populations on 
genetically differentiated than geographically distant 
ones. This was further supported by the fact that 
significantly less genetic differentiation and high 


average are more 


gene flow was observed between population pairs that 
were geographically further apart and vice versa. The 
presence of local barriers to the gene flow could be 
the reason for higher differentiation between nearer 
populations. The Vindhya mountain range might be 
acting as a barrier in the populations of Sagar and 
Umaria. In addition the vast forest cover could also 
affect the gene flow of urban malaria vector An. 
stephensi. The negative role of geographical distance 
in differentiation indicated the possibility of passive 
migration of An. stephensi which could be influenced 
by heavy human transport in these areas. 

The UPGMA dendrogram constructed according 
to Nei (1978) indicated that population of Gwalior, 
Sagar and Anuppur are closer genetically than 
Umaria (Fig. 6). 


Gwalior 
Sagar 
Anuppur 
Umaria 
eS 
2 


Fig. 6 UPGMA dendrogram showing genetic distance 
between populations of Anopheles stephensi 


The study revealed that the 
differentiation within eastern populations ( under 
EMCP) was higher as compared to that the northern 
populations (non-EMCP). Higher gene flow was 
seen within northern populations than the eastern. 


genetic 


However, significant gene flow was seen between 
EMCP and non-EMCP areas (Table 5), except the 
population of Umaria which seems to be isolated from 
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all other populations where less gene flow was 
observed. This could be ascribed to the stretches of 
Vindhya mountain range between these populations/ 
regions which might act as an effective barrier in 
gene flow. In addition, thick forest cover in and 
around the Umaria might also be preventing the gene 
flow. Further studies are needed to understand the 
role of Vindhya range to prevent the gene flow and 
transmission of malaria. In 1997, Enhanced Malaria 
Control Programme ( EMCP ) was launched by 
National Vector Borne Disease Control Programme 
(NVBDCP) with World Bank assistance to control 
tribal malaria. Despite World Bank initiative several 
focal outbreaks swept across the state between 2003 
and 2006. Though, the number of malaria infections 
and P. falciparum declined (77% and 69% , 
respectively ) initially as a result of intensive 
intervention measures in districts under EMCP, 
malaria cases started increasing (up to 50%) in 
areas which are not under EMCP and where malaria 
was never a problem. The number of deaths also 
showed an alarming increase (Singh, 2006). Our 
study revealed a high gene flow between the EMCP 
and non-EMCP areas of North Eastern Madhya 
Pradesh which coupled with the spread of insecticide 
resistance could be one of the main reasons for the 
cases, in spite of 
adequate control 
strategies. The development of resistance in An. 


increase in the malaria 


implementation of malaria 


culicifacies mosquito to pyrethroids has already been 
Pradesh. 


extensively used in malaria control programmes for 


reported in Madhya Pyrethroids are 
indoor residual spray and for impregnation of bed 
nets (Mishra et al., 2012). High migration rate in 
these areas could result in spread of insecticide 
resistance resulting in increase and stabilization of 
mosquito population. 
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WE: [HH RIKE Anopheles stephensi Ke WYN AR PA ABT A AE BE AY E BEEP, EN BE 12% PE BEI Bl HE S| 
fe, ASIF T EV BEA AB Madhya Pradesh Hb DX ZR JERR AJER IRETE il (EMCP ) XX AAR SR Le tll (IE EMCP ) Kr 
RIZICI TE EMCP XP RAT AAR ESE A EPO A E SEAR (AER RIAL DEW SE 
PRA or AEF RE [DEMAT “TB i HE BE AS Madhya Pradesh HY XX ZR IE AB AY 4 4+ EMCP X JE 
EMCP XR SE HUST RIEDE TAEK , VAT EE, [ERERKEN A aE PF EE AS 
HE, TERI RIAA UR A MIR EA RPE. WUE EMCP XW ARB RPEE( Pep = 0.0485, Rep =O. 1112) KJE EMCP X 
WAL ARE (Pop = 0.020, Rey = 0.0145) BA Bey IIE aM, HE EMCP KAJE EMCP K Z [i] WLS Bee AY EA Tit 
(12.90, 6.16, 5.06 12.38), RoW R BUET For, BLUME TT AE AE FA 8 2S TSE TE ee |. «Lave ) 8 
WIE eH , TE EMCP X Ade EMCP XLA EMCP Xd EMCP XIAP ER SSE A. SER CEE Rt 
EPER REFE EMCP XX AE EMCP KIERA RES EE 

Keil: WR. DB poids FA; WHEL; BEA; PERE 

FARIS: Q963 XRRR: A = EHH: 0454-6296 (2015 )01-0074-08 
































































































































































































































